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Classical  Multiscale  Analysis 

Wavelets:  Enormous  impact 

-  Theory 

-  Applications 

-  Many  success  stories 

Deep  understanding  of  the  fact  that  wavelets  are  not  good  for  all  purposes 
Consequent  constructions  of  new  systems  lying  beyond  wavelets 

Overview 

Other  multiscale  constructions 

Problems  classical  multiscale  ideas  do  not  address  effectively 
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Fourier  Analysis 

•  Diagonal  representation  of  shift-invariant  linear  transformations;  e.g. 
solution  operator  to  the  heat  equation 

dtu  =  a2  Am 

Fourier,  Theorie  Analytique  de  la  Chaleur,  1812. 

•  Truncated  Fourier  series  provide  very  good  approximations  of  smooth 
functions 


if  /  €  Ck. 


11/ -  Sn(f)\\L2  <C-n~k, 


Limitations  of  Fourier  Analysis 

Does  not  provide  any  sparse  decomposition  of  differential  equations  with 
variable  coefficients  (sinusoids  are  no  longer  eigenfunctions) 

Provides  poor  representations  of  discontinuous  objects  (Gibbs 
phenomenon) 


Wavelet  Analysis 

Almost  eigenfunctions  of  differential  operators 

(■ Lf)(x )  =  a(x)dxf(x) 


Sparse  representations  of  piecewise  smooth  functions 


Wavelets  and  Piecewise  Smooth  Objects 

1 -dimensional  example: 

9(t)  =  1  {t>Me_(t_to)2- 

Fourier  series: 


g(t)  =  £  Ckeikt 

Fourier  coefficients  have  slow  decay: 

|c|(„)  >  c  •  l/n. 

Wavelet  series: 

g(t)  =  £6>a  xpx(t) 

Wavelet  coefficients  have  fast  decay: 
l0l(n)  <  c  •  (l/n)r  foranyr>0. 


as  if  the  object  were  non-singular 
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’Wavelets  and  Operators’ 

Calderon-Zygmund  Operator 

(■ Kf)(x)  =  J  K(x,y)f(y)dy 

•  K  singular  along  the  diagonal 

•  K  smooth  away  from  the  diagonal 

•  Sparse  representation  in  wavelet  bases  (Calderon,  Meyer,  etc.) 

•  Applications  to  numerical  analysis  and  scientific  computing  (Beylkin, 
Coifman,  Rokhlin) 


Our  Viewpoint 

Sparse  representations  of  point-singularities 
Sparse  representation  of  certain  matrices 
Simultaneously 

Applications 

-  Approximation  theory 

-  Data  compression 

-  Statistical  estimation 

-  Scientific  computing 

More  importantly:  new  mathematical  architecture  where  information  is 
organized  by  scale  and  location 


New  Challenges 

Intermittency  in  higher-dimensions 
Evolution  problems 

CHA  has  not  addressed  these  problems. 


Agenda 

Limitations  of  existing  image  representations 
Curvelets:  geometry  and  tilings  in  Phase-Space 
Representation  of  functions,  signals 
Representation  of  operators,  matrices 


Three  Anomalies 

Inefficiency  of  Existing  Image  Representations 


Limitations  of  Existing  Pyramid  Schemes 
Limitations  of  Existing  Scaling  Concepts 
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I:  Inefficient  Image  Representations 

Edge  Model:  Object  /( a?i,  x2 )  with  discontinuity  along  generic  C2  smooth 
curve;  smooth  elsewhere. 

Fourier  is  awful 

Best  m-terrn  trigonometric  approximation  fm 


ll/-/m.|l2Xm  1/2>  rn  — ►  OO 

Wavelets  are  bad 

Best  m-term  approximation  by  wavelets: 


11/  -  fm 


—  1 


2  X  m 


5 


rn 


oc 
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Optimal  Behavior 

There  is  a  ‘dictionary’  of  ‘atoms’  with  best  m-term  approximant  fm 


11/  ~  fm 


m 


oo 


•  No  basis  can  do  better  than  this. 

•  No  depth-search  limited  dictionary  can  do  better. 

•  No  pre-existing  basis  does  anything  near  this  well. 


(a)  Wavelets 
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(b)  Triangulations 


II.  Limitations  of  Existing  Pyramids 

Canonical  Pyramid  Ideas  (1980-present) 

•  Laplacian  Pyramid  (Adelson/Burt) 

•  Orthonormal  Wavelet  Pyramid  (Mallat/Meyer) 

•  Steerable  Pyramid  (Adelson/Heeger/Simoncelli) 

•  Multiwavelets  (Alpert/Beylkin/Coifman/Rokhlin) 

Shared  features 

•  Elements  at  dyadic  scales/locations 

•  FIXED  Number  of  elements  at  each  scale/location 
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Wavelet  Pyramid 


Limitations  of  Existing  Scaling  Concepts 

Traditional  Scaling 

fa(xi,x2)  =  f(ax1,ax2),  a  >  0. 


Curves  exhibit  different  kinds  of  scaling 

•  Anisotropic 

•  Locally  Adaptive 

If  f  (x x2 )  —  then 


fa(x i,x2)  =  f(a  •  x1,a2x2) 


In  Harmonic  Analysis  called  Parabolic  Scaling. 
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Figure  1 :  Curves  are  Invariant  under  Anisotropic  Scaling 


Curve  lets 


C.  and  Guo,  2002. 

New  multiscale  pyramid: 

•  Multiscale 

•  Multi-orientations 

•  Parabolic  (anisotropy)  scaling 

width  ~  length 


Earlier  construction,  C.  and  Donoho  (2000) 
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Philosophy  (Slightly  Inaccurate) 

•  Start  with  a  waveform  ip(x )  =  cp(x i,  x2). 

-  oscillatory  in  a?i 

-  lowpass  in  x2 

•  Parabolic  rescaling 


Dj\ip(Djx)  =  23j’/4i^(2J'a:i,  2^2x2), 


0 

2^'/2 


,  J  >0 


•  Rotation  (scale  dependent) 

23^4^(DjR0jtx),  Oji  =  2t t  •  £2~ U/2J 

•  Translation  (oriented  Cartesian  grid  with  spacing  2~j  x  2~j/2)\ 


2  V^ipiDjRo  x  -  fc),  fee  Z2 
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Curvelet  Construction,  I 

Decomposition  in  space  and  angle 

•  ‘Orthonormal’  partition  of  unity  (scale):  smooth  pair  of  windows  (w0,  wi) 

oo 

wo(r)  +  =  1?  r  >  0. 

j=o 


•  ‘Orthonormal’  partition  of  unity  (angle): 

oo 

^  v2(t  —  £)  =  1,  t  G  M. 

£=  —  oo 

•  For  each  pair  J  =  (j,  £),  define  the  two-dimensional  window 

Wj(£)  =  «;i(2-^|i|)  •  v(2VW(0  -  OjJ),  0jfe  =  tt  •  £  ■  2~V/*i. 


Note 


Wj  is  an  orthonormal  partition  of  unity 


53  Wj(£)  =  1- 

J 

Interpretation 


Divide  frequency  domain  into  annuli  |a?|  e  [2j ,  2j+1) 
Subdivide  Annuli  into  wedges.  Split  every  other  scale. 
Oriented  local  cosines  on  wedges 
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Partition:  Frequency-side  Picture 


Curvelet  Construction,  II 

Local  bases 

•  Fix  scale/angle  pair  (j,£) 

suppW^  C  Dd,e  =  Re'D^o 

where  Dj>0  is  a  rectangle  of  sidelength  2w2j  x  2Tz2j/2,  say. 

•  Rectangular  grid  A j  of  step-size  2~j  x  2~j I2 

A j  =  {k:k  =  (k12-j,k22~j/2),k1,k2  E  1} 


Local  Fourier  basis  of  L2(Dj,0) 


ej,  0,fe(C)  — 


27t23j/4 


k  €  A7- 


and  likewise  local  Fourier  basis  of  L2(Dj >e) 


ej,£,k  —  ej,o,k(Rej^O 
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•  Curvelets  (frequency-domain  definition) 


Properties 

Reproducing  formula 

/  =  y^(/?  ifj,k)‘pj,k 

J,k 

Why?  Let  g  be  the  RHS  and  apply  Parseval 

a  =  ^{f,vj,k)vj,k 

J,k 

J,k 

=  J2WJ^2(fWJ^J,k)ej,k 
J  k 

=  Y,  fwJwJ  =  fYw3  = 

j  j 

Conclusion  g  =  f  and,  therefore,  f  =  g. 

Parseval  relation  (similar  argument) 


/II2  =  Y  K/»  vj,k) 
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Curvelet:  Space-side  Viewpoint 

In  the  frequency  domain 

2~3o/4 

<Pj,0,k(Z)  =  - Wj, 0«)ei<fc»«,  k  e  Aj 

Z7T 

In  the  spatial  domain 

<Pj,o,k(x)  =  23j/4cpj(x  -  k),  Wj, o  =  2 it cp j 

and  more  generally 

lPj,e,k(x)  =  2 33/4ipj(Rgji(x  -  Rg'k)), 

All  curvelets  at  a  given  scale  are  obtained  by  translating  and  rotating  a  single 
‘mother  curvelet.’ 
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Further  Properties 


•  Tight  frame 

f  =  'y  '  (/ ’  V’j/jfc)  ^Pj,e.,k 

3  ,£,fc 


53</>  ^,fe)2 


•  Geometric  Pyramid  structure 

-  Dyadic  scale 

-  Dyadic  location 

-  Direction 

•  New  tiling  of  phase  space 


Needle-shaped 
Scaling  laws 

-  length  ~  2_J/2 

-  width  ~  2~ * 


width  length 


-  # Directions  =  2U/2J 

-  Doubles  angular  resolution  at  every  other  scale 
Unprecedented  combination 
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Second  Dyadic  Decomposition 

•  Fefferman  (70’s),  boundedness  of  Riesz  spherical  means 

•  Stein,  and  Seeger,  Sogge  and  Stein  (90’s),  Lp  boundedness  of  Fourier 
Integral  Operators. 

•  Smith  (90’s),  atomic  decomposition  of  Fourier  Integral  Operators. 
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Digital  Curvelets 


Source:  DCTvUSFFT  (Digital  Curvelet  Transform  via  USFFT’s),  C.  and  Donoho  (2004). 
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Digital  Curvelets:  Frequency  Localization 


50 

100 

150 

200 


i 


50  100  150  200  250  300  350  400  450  500 


Time  domain 


50 

100 

150 


•iso  If 

50  100  150  200  250  300  350  400  450  500 

Frequency  domain  (modulus) 


Source:  DCTvUSFFT  (Digital  Curvelet  Transform  via  USFFT’s),  C.  and  Donoho  (2004). 
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Digital  Curvelets:  :  Frequency  Localization 


Time  domain  Frequency  domain  (real  part) 


Source:  DCTvUSFFT  (Digital  Curvelet  Transform  via  USFFT’s),  C.  and  Donoho  (2004). 
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Curvelets  and  Edges 
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Optimality 

•  Suppose  /  is  smooth  except  for  discontinuity  on  C2  curve 

•  Curvelet  m-term  approximations,  naive  thresholding 


11/  -  / 


curve 

rn 


2  <  C  m  2(logm)3 


Near-optimal  rate  of  m-term  approximation  (wavelets  ~  m_1). 


Idea  of  the  Proof 


Decomposition  of  a  Subband 
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s 

Frequency  2 


/ 
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Edge  Fragments 

Parabolic  edge  Square  root  of  the  truncated  absolute  value  of  the  FT 


50 

100 

150 

200 

250 

300 

350 

400 

450 

500 

50  100  150  200  250  300  350  400  450  500  50  100  150  200  250  300  350  400  450  500 


Edge  fragment 


FT  of  an  edge  fragment  (abs) 
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Curvelets  and  Warpings 

C2  change  of  coordinates  preserves  sparsity  (C.  2002). 


•  Let  :  R2  — ►  R2  be  a  one  to  one  C2  function  such  that  \\JV 
bounded  away  from  zero  and  infinity. 


oo 


is 


•  Curvelet  expansion 


f  =  Y1  M/)7m>  M/)  =  </»  7m) 

M 


•  Likewise, 

f  °<P  =  °  ¥>)7m 

M 

The  coefficient  sequences  of  /  and  /  o  ip  are  equally  sparse. 
Theorem  1  Then ,  for  each  p  >  2/3,  we  have 


I|0(/O¥>)|UP  <Cp.|W)H/p. 
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Curvelets  are  nearly  invariant  through  a  smooth  change  of  coordinates 
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Curvelets  and  Curved  Singularities 

•  /  smooth  except  along  a  C2  curve 

•  fn,  n-term  approximation  obtained  by  naive  thresholding 

11/  -  /n|li2  <  C  •  (log  n)3  •  n~ 2 


•  Why? 

1 .  True  for  a  straight  edge 

2.  Deformation  preserves 
sparsity 

•  Optimal 


<i> 


arbitrary  singularity 


straight  singularity 
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Curvelets  and  Hyperbolic  Differential  Equations 


Representation  of  Evolution  Operators 

Wave  Equation 

d2u  =  c2(x)Au,  x  G  R2,  R3,  ... 


Evolution  operator  S(t ;  x,  y )  is  dense! 

Modern  viewpoint 

-  Basis  (pn  of  L2. 

-  Representation  of  S :  A(n,n')  =  (v?n,  Syjn/) 


u  =  <Stio 


F 

d(u0) 


0{u)  =  S0(no) 


Find  a  representation  in  which  S  is  sparse 
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Range  of  Influence 


The  solution  operator  S  is  a  ‘dense’  integral. 
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Potential  for  Sparsity 

If  the  matrix  A  is  sparse,  potential  for 

•  fast  multiplication 

•  fast  inversion 


Example:  convolutions  and  Fourier  transforms 
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Current  Multiscale  Thinking 

•  Traditional  ideas: 

-  Multigrids 

-  Wavelets 

-  Adaptive  FEM’s 

-  etc. 

•  Traditional  multiscale  ideas  are  ill-adapted  to  wave  problems: 

1 .  they  fail  to  sparsity  oscillatory  integrals  like  the  solution  operator 

2.  they  fail  to  provide  a  sparse  representation  of  oscillatory  signals  which 
are  the  solutions  of  those  equations. 


Peek  at  the  Results 


Sparse  representations  of  hyperbolic  symmetric  systems 
Connections  with  geometric  optics 
Importance  of  parabolic  scaling 


Symmetric  Systems  of  Differential  Equations 

dtu+Y^  A j ( t/* ^ |  B(x)u  —  3? J  —  i  ®  G  R 

i 

u  is  m-dimensional  and  the  Ai  s  are  symmetric 
Examples 

•  Maxwell 

•  Acoustic  waves 

•  Linear  elasticity 
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Example:  Acoustic  waves 

System  of  hyperbolic  equations:  v  =  ( u,p ) 


•  Dispersion  matrix 


dtu  +  V  (c(x)p) 

— 

0 

dtp  +  c(x)  V  •  u 

— 

0 

( 0 

0 

ki\ 

Y^Aj(x)kj  =  c(x) 

0 

0 

k2 

3 

\fel 

fe2 

0/ 

•  Eigenvalues  (A„):  A±  =  ±c(a:),  A0  =  0. 

•  Eigenvectors  (Ru) 


R0(k)  = 


R±(k)  = 


1 

V5 


Geometric  Optics:  Lax,  1957 


High-frequency  wave-propagation  approximation 


V 


(x,t)  =  f aou(x,t )  + 

is  ' 


t)  t) 


+ 


+ 


u; 


•  Plug  into  wave  equation 

-  Eikonal  equations 

dt&v  +  \„(x,  =  0. 

\v(x,  k)  are  the  eigenvalues  of  the  dispersion  matrix  V .  Aj(x)k 

-  And  ’transport’  equations  for  amplitudes 

Turns  a  linear  equation  into  a  nonlinear  evolution  equation! 
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Hamiltonian  Flows 

•  Dispersion  matrix 

y  ]  Aj(x)kj 

3 

•  Eigenvalues  Xu(x,  k),  eigenvectors  Ru{x,  k ) 

•  Hamiltonian  flows  (in  general,  m  of  them)  in  phase-space 

(x{t)  =Vk\u(x,k),  cc(0)  =  x0, 

k(t)  =  —Vx\v(x,k),  fc(0)  =  fe0. 

•  Eikonal  equations  from  geometric  optics 

dt&L,  +  \„(x,  Va;#)  =  0. 

$  is  constant  along  the  Hamiltonian  flow  $(t,  x(t))  =  Cste.  Problem: 
valid  before  caustics,  i.e.  before  <f>  becomes  multi-valued  (ray-crosses). 
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Hyper-curvelets 


•  Frequency  definition 

=  Ejycpll(k). 

•  Hyper-curvelets  build-up  a  (vector-valued)  tight-frame,  namely 


u 


2 

L2 


—  I 


2 

5 


and 


•  Other  possibilities: 


^  1//X.  (fc)  =  R„(k)v„(k), 


R„(x,k)ip^(k)ei{k’x)  dk. 


Approximation 

The  action  of  the  wave  propagator  on  a  curvelet  is  well-  approximated  by 
rigid  motion. 

•  Initial  data:  x ) 

•  Approximate  solution 

0 t ,  fe)  =  R„(k)<p^t)(k) 


where 

W(t)(«)  =  W  (C^/x(*)(a5  -  a5M(t))  +  ®M) 

-  is  a  rotation 

-  Xfj,(t)  is  a  translation 
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Main  Result 

Curvelet  representation  of  the  propagator  S(t)\ 

A(v,p-,v',p')  =  S(t)&u>^) 

Theorem  2  (C.  and  Demanet)  Suppose  the  coefficients  of  a  general 
hyperbolic  system  are  smooth,  i.e.  C°°. 

•  The  matrix  is  sparse.  Suppose  a  is  either  a  row  or  a  column  of  A,  and  let 
|a|(Tl)  be  the  n-largest  entry  of  the  sequence  |a|,  then  for  each  r  >  0, 
|a|(n)  obeys 

1^1(71.)  ^  CrTi  . 

•  The  matrix  is  well-organized.  There  is  a  natural  distance  over  the  curvelet 
indices  such  that  for  each  M, 

\a(is,  p;  i/,  p')\  <  CM,t(vi  v')  ^(1  +  d(p,  pu»(t))~M 


The  constant  is  growing  at  most  like  C\eC2t,  and  for  v  ^  v' ,  Ci  <<  1. 
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Sketch  of  the  curvelet  representation  of  the  wave  propagator 
A(is',  n'\  v,  n)  =  S (t)<& . 
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Our  Claim 

Curvelets  provide  an  optimally  sparse  representation  of  wave  propagators. 

•  Fourier  matrix  is  dense 

•  FEM  matrix  is  dense 

•  Wavelet  matrix  is  dense 

Second-Order  Scalar  Equations 


d^u  —  a,ij(x)didjU  =  0 

Sparsity  comes  for  free  (in  the  scalar  curvelet  system) 
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Why  Does  This  Work? 

“The  purpose  of  calculations  is  insight,  not  numbers...” 

(adapted  from  R.  W.  Hamming) 


•  Wave  operator  (constant  coefficients) 

d2u  =  c2Au,  u(x,  0)  =  u0(x),  dtu(x,  0)  =  «i(a;). 


•  Fourier  transform 

u(t,  k )  =  J  u(t,  x)e~ tk'x  dx. 

•  Solution 


.„n  ,  I,  \  ,  Sin(c|fe|t)  . 

u(t,k)  =  cos(c\k\t)Uo(k) -\ - - — - - Ui(k) 


_  g=bic|fc|t 


|fc| 

.zbicl  felt 


u0(k)  + 


\k\ 


ui(k ) 


■u(t,  k) 


Set  tr  =  ct 


ei\k\t'  _  eiklt'+S(k)t _  |fc|-fci. 

Because  of  parabolic  scaling 

S(k)  =  0(1) 


Frequency  modulation  is  nearly  linear 


ei\k\f  _  ^ k ^ .  eiklt'' 


h  smooth  and  non-oscillatory 
Space-side  picture: 

1 .  modulation  corresponds  to  a  displace¬ 
ment  of  ±ct  in  the  codirection. 

2.  multiplication  by  h\  gentle  convolution 
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cos(|k|  t')  ~  cos^t') 


Displacement  in  the 
codirection 
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Displacement  in  the 
codirection 


Warpings 

Variable  coefficients:  utt  =  c2(x)Au. 

At  small  scales,  the  velocity  field  varies  varies  little  over  the  support  of 
curvelet. 

Model  the  effect  of  variable  velocity  field  as  a  warping  g\  <p^(g(x)). 
Warpings  do  not  distort  the  geometry  of  a  curvelet. 
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Importance  of  Parabolic  Scaling 

Consider  arbitrary  scaling  (anisotropy  incresases  as  a  decreases) 

width  ~  2~^ ,  length  ~  2~-’°l,  0  <  a.  <  1. 

-  ridgelets  a  =  0  (very  anisotropic), 

-  curvelets  a  =  1/2  (parabolic  anisotropy), 

-  wavelets  a  =  1  (roughly  isotropic). 

For  wave-like  behavior,  need 

0 

width  <  length 

For  particle-like  behavior,  need 

0  o 

width  >  length 

For  both  (simultaneously),  need 

0  r\ 

width  length 


Only  working  scaling:  a  =  1/2 


Examples  I 


t  =  0 


A 
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Examples  II 


Examples  III 


t  =  o 


t  =  T 
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Architecture  of  the  Argument 

1 .  Decoupling  into  polarized  components:  Decompose  the  wavefield  u(t,  x ) 
into  m  one-way  components  fu(t ,  x) 

rn 

u(t,x)  —  'y  ^  Rvfv{t,  a;), 

i/=i 

Ru\  pseudo-differential  operator  (independent  of  time),  maps  scalar  fields 
into  m-dirnensional  vector  fields. 

rn 

S(t )  =  ^2  R^e~'ltAl' L„  +  negligible, 

U=l 

Lu:  pseudo-differential  operator,  maps  m-dirnensional  vector  fields  into 
scalar  fields. 
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Curvelet  Splitting 


Time  t  =  0 


Time  t 
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s(t)  ~  -M*)*  =  R^~itKvLv 

2.  Fourier  Integral  Operator  parametrix.  We  then  approximate  for  small  times 
t  >  0  each  e~itKv ,  v  =  1, . . . ,  m,  by  an  oscillatory  integral  or  Fourier 
Integral  Operator  (FIO)  F„(t) 

Fv(t)f(x)  =  J  ez*,'(t,x’k'> cr„(t,  x,  k)  f(k)  dk 

Key  issue :  decoupling  is  crucial  because  this  what  makes  it  work  for  large 
times  (not  an  automatic  consequence  of  Lax’s  construction) 

S„(nt)  =  S„{t)n 

3.  Sparsity  of  FIO’s.  General  FIO’s  F(t)  are  sparse  and  well-structured  when 

represented  in  tight  frames  of  (scalar)  curvelets  — a  result  of 

independent  interest. 


Potential  for  Scientific  Computing 

Transform  this  theoretical  insight  into  effective  algorithms 

Sources  of  inspiration 

-  Rokhlin 

-  Engquist  &  Osher 

-  Others 
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Our  Viewpoint 


ut 


e 


-pt 


u  o 


—pt 


Uq 


->  Ut 


#0 


A(t) 


+  e, 


For  any  t,  A(t )  is  sparse 


Background:  Fast  and  accurate  Digital  Curvelet  Transform  is  available  (with 
D.  Donoho). 


Grid  size  N 


Preliminary  Work 


Accuracy  e 
Complexity 

0(N1+s-e~d),  any  <5  >  0 


Conjecture 


0(N  log  N  •  e~d) 


Arbitrary  initial  data 
Works  in  any  dimension 
General  procedure 
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Propagating  curvelets  is  not  geometric  optics! 


Wave  Front,  initially 


Summary 

New  geometric  multiscale  ideas 
Key  insight:  geometry  of  Phase-Space 
New  mathematical  architecture 
Addresses  new  range  of  problems  effectively 
Promising  potential 


Numerical  Experiments,  I 
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(a)  Noisy  Phantom 
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(b)  Curvelets 


(c)  Curvelets  and  TV 


(d)  Curvelets  and  TV:  Residuals 
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(e)  Wavelets 


(f)  Curvelets 
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(I)  Noisy 
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(m)  Curvelets  &  TV:  Residuals 


(n)  Noisy  Scanline 
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(o)  True  Scanline  and  Curvelets  and  TV  Re¬ 
construction 


(a)  Noisy  Scanline 


1.5 


(c)  Curvelets  and  TV 


(b)  Curvelets 


(d)  True  Scanline  and  Curvelets  and  TV  Re¬ 
construction 


(p)  Original 


(q)  Noisy 
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(t)  Noisy  Detail 


(u)  Curvelets 


(v)  Curvelets  and  TV 


Numerical  Experiments,  II:  Work  of  Jean-Luc  Starck  (CEA) 


a)  Simulated  image  (Gaussians+lines)  b)  Simulated  image  +  noise 


c)  A  trous  algorithm 


d)  Curvelet  transform 


e)  coaddition  c+d 


f)  residual  =  e-b 


Separation  of  Texture  from 
Piecewise  Smooth  Content 


The  separation  task:  decomposition,  of  an  image 
into  a  texture  and  a  natural  (piecewise  smooth; 
scene  part. 


